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We present the first successful application of the method of Matched Expansions for the 
calculation of the self-force on a point particle in a curved spacetime. We investigate the 
case of a scalar charge in the Nariai spacetime, which serves as a toy model for a point 
mass moving in the Schwarzschild black hole background. We discuss the singularity 
structure of the Green function beyond the normal neighbourhood and the interesting 
I , effect of caustics on null wave propagation. 
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1. Introduction 

Extreme Mass Ratio Inspirals, such as a stellar- mass compact object orbiting around 
a supermassive black hole, can be understood within perturbation theory of General 
ly^ . Relativity as the smaller body's motion deviating from the background geodesic due 

to a self -force. 1 The calculation of the self- force is necessary in order to model the 
■ orbital motion and obtain accurate gravitational- wave templates. The expression for 

, the self-force in a general spacetime is formally similar for the cases of an orbiting 

point scalar charge, electrical charge and a point mass. In all these cases, the self- 
force is given by some local terms which are relatively easy to calculate, plus a 'tail' 
term which, in the case of a scalar charge q following a worldline z(t), is given by 



q 2 f V^G ret (z(T),z(T'))dT' (1) 

J — CO 

where G re t(x, x') is the 'retarded' Green function of the scalar field wave equation. 
As the 'tail' term depends on the entire past history of the particle's worldline, it 
is usually quite difficult to calculate and one has to do so numerically. 

2. Method of Matched Expansions 

We choose to use the method of Matched Expansions 2 to calculate (pQ). In this 
method, the 'tail' integral is separated into two regimes: 

(a) In the quasilocal regime, r' G (r — At, r) for a certain 'small value' At, the 
Green function is calculated using the Hadamard form: 



G 



ret 



(x, x') = 6L (x, x') {A 1 / 2 (x, x') 5{a)-V (x, x') 6 (-<r)} (2) 
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where a is half the square of the geodesic distance between the spacetime points x 
and x 1 ', the van Vleck determinant A(x,x') and V(x,x') are two regular biscalars 
and 0- (x,x f ) guarantees causality. Given the structure of (j2j) it is clear that only 
V(x, x') will contribute to the 'tail' integral in the quasilocal regime. The Hadamard 
form is only valid for points 'nearby' - within a normal neighbourhood. 

(b) In the distant past regime, r' G (— oo,r — At), the Green function can be 
calculated using a Fourier mode decomposition and deforming the frequency-integral 
into a contour in the complex frequency-plane. 

3. Singularity structure of the Green function 

We carried out 3 the first succesful application of the method of Matched Expan- 
sions. With the future view of applying it to the case of a point mass moving 
in an astrophysically-relevant spacetime, we used a toy model where calculations 
are easier: a scalar charge moving in the Nariai spacetime (cartesian product of 
2-dimensional de Sitter spacetime and the 2-sphere), in which the radial potential 
of Schwarzschild is replaced by the Poschl- Teller potential. 

In the quasilocal regime, some powerful methods have been developed 4 ' 5 for the 
calculation of V(x,x') up to, and beyond, the normal neighbourhood. 

In the Nariai spacetime, the 'retarded' Green function in the 'distant past' is 
completely given by the sum over its poles - the quasinormal modes. A study of this 
series beautifully uncovered - FigJT]- the singularity structure of the Green function: 
its singularities occur at the times it takes a null geodesic to travel between x and 
x' (i.e., when a = 0); furthermore, the singularity structure follows a 4- fold cycle: 



An alternative, heuristic way for understanding how this structure arises is by 
use of the Feynman Green function, which has Hadamard form 



where W(x, x') is another regular bitensor and e > 0. The 'retarded' Green function 
can be obtained as G ret (x, x r ) = 26L(x, x') lim e ^o Re (Gf(x, x')). Although strictly 
speaking the Hadamard form is only valid within the normal neighbourhood, if we 
tentatively use it beyond we can obtain the van Vleck determinant from a transport 
equation past caustic points, where it diverges. One then finds that A 1 / 2 (x,x / ) 
picks up a factor '— V every time a null geodesic crosses a caustic. From the property 
lim e ^ + i/(o~-\-ie) = 7r5(a) + i/o~, one can see how the 4- fold cycle arises. This factor 
seems to be typical of the S 2 topology and, although it is known in other fields of 
Science, in General Relativity it has only been previously observed by Prof. Ori, 6 
to the best of our knowledge. 

We were able to successfully match the quasilocal and 'distant past' contribu- 
tions to the 'retarded' Green function - Figj2^a). The calculation of the 'retarded' 
Green function over the whole past of the particle's worldline enabled us to calculate 
the self- force for a scalar charge in Nariai - Figj2^b). 



8(a), I /(va), -8(a), -I /(va), 5(a),.... 
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Fig. 1. (a) Null geodesies orbiting the unstable photon orbit of the background, (b) 'Retarded' 
Green function as a function of the proper time of the static scalar charge in Nariai. There is 
a 1-to-l correspondence between the geodesies and the times when the Green function becomes 
singular, with a four- fold structure: S(a), — , —5(a), — — , S(a), . . . 



Matching Quasi-local and Non-local Green Functions (f=l/8) 



Radial Self-Force on a Static Particle in Nariai Spacetime 



- Quasilocal (60th order) 








• Smoothed Sum (n=0) 






0.0004 


Watson Transform (n=0) 














IK 








m 0.0003 








M 






| 0.0002 


















1 






0.0001 




0.0000 




Fig. 2. (a) Matching of the 'quasilocal' and 'distant past' contributions as a function of time, 
(b) self- force for a static scalar charge in the Nariai spacetime as a function of the radius. 
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